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1.  Contrast  Otyectives 

The  objective  of  the  contract  was  to  extend  the  work  of  the  earlier  feasibility 
study[1]  on  the  control  of  the  operating  frequency  of  milliinetre-wave  IMPATT 
oscillators  using  an  optical  signal.  Specific  objectives  were  to  maximise  the 
frequency  tuning  obtainable  with  conventional  IMPATT  structures,  to  fabricate  a 
special  IMPATT  structure  designed  for  optical  control  and  to  examine  tuning 
speed  and  ncdse  performance. 

2.  Summary  of  Principal  Results 

The  analytic  theory  of  optical  tuning  in  IMPATT  oscillators  developed  in  the 
initial  phase  of  the  work  [1]  has  been  extended  and  generalised.  Accuracy  of 
the  theory  in  predicting  tuning  at  the  higher  oscillator  voltage  swings  has  been 
greatly  improved  by  reformulating  the  Bessel  function  series  used  to  express 
avalanche  gain.  A  model  for  optical  injection  locking,  where  the  oscillator 
output  is  locked  to  a  microwave  intensity  modulated  optical  signal,  has  also  been 
developed. 

The  computer  model  for  optical  control  has  been  modified  to  study  the  tuning 
behaviour  of  Ku-Band  IMPATTs,  as  used  in  the  recent  experimental  work.  Work 
on  optical  tuning  speed  studies  has  also  been  extended. 

The  experimental  programme  has  produced  a  number  of  developments.  A  delay 
line  discriminator  has  been  used  to  measure  the  tuning  speed  in  response  to  a 
step  change  in  illumination  level-  It  has  been  found  that  a  frequency  shift  of 
1.4MHz  at  W-Band  could  be  achieved  in  less  than  20ns  and  a  frequency  shift  of 
5.6MHz  at  Ka-Band  in  less  than  10ns.  However  both  of  these  results  are  thought 
to  be  due  to  limitations  in  the  measurement  system  rather  than  in  the  intrinsic 
response  time.  Nonetheless  these  results  are  valuable  in  that  they  confirm  that 
the  tuning  response  is  not  associated  with  trapping  or  other  slow  effects. 

A  special  Ka-Band  IMPATT  structure  was  fabricated,  using  a  modified  integral 

♦ 

heat  sink  process  to  produce  a  40  lira  diameter  hole  beneath  the  P  region  of  the 

i . 


diode  through  which  the  diode  could  be  illuminated.  Comparative  teats  showed 
that  the  tuning  slope  was  improved  with  through  heat-sink  illumination  by  almost 
an  order  of  magnitude. 


Finally  modifications  to  edge  illumination  techniques  for  IMPATT  diodes  have  been 
investigated  including  the  use  of  quartz  stand-off  rather  than  quartz  ring 
mounting.  Also  a  technique  for  reducing  the  diameter  of  optical  fibres  by  wet 
chemical  etching  has  been  developed,  which  has  applications  in  a  number  of 
illumination  configurations. 


3-  Analytic  Theory 

When  the  electric  field  in  an  IMPATT  diode  is  above  the  breakdown  value  charge 
carriers,  whether  optically  generated  or  not,  experience  avalanche  multiplication. 
This  avalanche  multiplication  was  represented  in  the  earlier  theory  [1]  by 
truncated  Bessel  function  series  of  argument. 


X  = 


2(ra+1)V 


_L 


•(1) 


un  V 
a  0 


where  m  is  an  ionisation  non-linearity  constant  (ms  6  for  silicon), is  the  peak 

value  of  the  modulated  component  of  the  voltage  across  the  IMPATT, V  is  the 

b 

breakdown  voltage, u>  the  operating  frequency  and  Tg  the  avalanche  zone  transit 
time,  and  where  It  is  assumed  that  the  device  remains  punched  through  over  the 
whole  cycle. 


As  the  vcdtage  swing  across  the  device  increases  higher  order  BesseL  function 
coefficients  become  significant  and  the  truncated  series  give  a  poorer 
representation  of  the  effect  of  avalanche  multiplication.  The  extended  analytic 
theory  uses  more  general  expressions  which  are  capable  of  good  accuracy  over  a 
wide  range  of  voltage  swings. 

At  the  same  time  the  opportunity  was  taken  to  extend  the  theory  to  treat 
optical  injection  locking. 
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In  optical  injection  locking  the  optical  control  signal  is  intensity  modulated  at  a 
frequency  close  to  the  free-running  frequency  of  the  IMPATT  oscillator(2]. 
Although  the  technique  has  yet  to  be  demonstrated  at  millimetre- wave 
frequencies  due  to  the  difficulty  of  intensity  modulating  the  control  signal  at 
such  high  rates,  developments  in  optical  heterodyne  technctlogy[3]  suggest  that 
the  technique  could  offer  novel  systems  possibilities  in  the  future. 

The  derivation  of  the  extended  theory  is  given  in  Appendix  A;  only  the  most 
important  results  will  be  discussed  here.  Considering  first  optical  tuning:  the 
change  in  oscillation  frequency  for  a  small  change  in  optically  generated  current, 
is  given  by: 
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where  Q  is  the  oscillator  Q,  I  the  oscillator  bias  current  and  F  {  X }  is  an 

dc  a 

avalanche  gain  term.  Comparison  with  equation  4.26  of  [1]  shows  that  the  only 
difference  lies  in  the  definition  of  the  avalanche  gain  parameter  F  f  x)  • 


Figure  la  shows  the  dependence  of  F  I  xl  on  X  as  an  increasing  number  of 

a 

terms  in  the  Bessel  function  series  defining  F  ( X}  (see  equation  10  of  Appendix 

a 

A)  are  taken  into  account.  To  appreciate  the  practical  significance  of  these 
curves  X  must  be  related  to  the  oscillator  voltage  swing  through  equation  (1). 
For  a  typical  uniformly  doped  device  Ta  ~  0.2  t  /  which  with  m  =  6 
gives 


X 


(3) 


The  earlier  theory  corresponds  to  the  curve  n  ^  1:  this  becomes  seriously  in 
error  for  X  >  3.0,  corresponding  to  >  14  <  vcQtage  modulation.  Curves  for  n  ^  2 
or  3  give  satisfactory  accuracy  for  up  to  25 1  voltage  modulation,  the  highest 


value  likely  to  be  used  in  uniformly  doped  structures. 


For  optical  injection  locking  the  full  locking  range  is  given  by 
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where  I  lx^  are  modified  Bessel  functions  and  F  f  X  lis  an  avalanche  gain 
n  b 

parameter.  The  optical  injection  locking  range  is  proportional  to  the  peak 

optically  injected  locking  current  I  .  Figure  1b  shows  the  dependence  of  F  { X  } 

.  Sj  b 


on  X.  From  equation  (3)  F  I  X 


250  for  20?  voltage  modulation;  this 


avalanche  gain  renders  IMPATT  diodes  particularly  suitable  for  optical  control. 


As  a  numerical  example,  consider  a  Ka-Band  oscillator  having  to  :  2"  x  3-8  x 
10  -1  o 

10  r  s  ,  Q  =  100  and  I  =  100mA. 

dc 


With  20?  voltage  modulation  and  a  peak  optically  injected  locking  current  of 
100  pA  the  predicted  flail  locking  range  would  be  540MHz,  a  practically  useful 
value. 


4.  Computer  Modelling  Studies 

Because  of  the  decision  to  broaden  the  work  to  include  Ka-Band  IMPATTs  the 
parameters  of  the  model  had  to  be  changed.  However  the  basic  model  described 
in  section  5-1  of  [1]  was  retained. 

4.1.  Device  parameters 

The  Ka-Band  devices  used  in  the  experlm  ental  work  were  silicon, 
♦  + 

uniformly  doped  P  NN  structures.  The  material  parameters  of  [1]  were 
retained,  with  an  assumed  Junction  temperature  of  450K.  A 
complementary  error  function  model  for  the  doping  profile  gave  a 


satisfactory  match  to  the  measured  values  using  the  parameters  listed 
below  in  equation  (5-7)  of  [1]: 
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Figure  2  shows  the  modelled  profile  and  Figure  3  the  measured  profile; 
agreement  is  seen  to  be  satisfactory. 

4.2.  DCP  Results 

-2 

A  bias  current  density  of  4,000  A  cm  was  selected  for  the  modelling 
work,  corresponding  to  a  bias  current  of  112mA  in  the  experimental 
devices. 

Figure  4  shows  the  electric  field  and  electron  current  density  profiles 

for  the  device.  The  diode  is  seen  to  be  punched-through,  thus 

minimising  losses  in  undepleted  material.  Applying  the  20%-80%  criterion 

to  the  electron  current  density  curve  suggests  an  avalanche  zone  length, 

L  =  0.23wn. 
a 

The  predicted  bias  voltage  was  29. 5 V  with  the  diode  not  oscillating,  a 
value  in  satisfactory  agreement  with  measured  values  for  the 
experimental  oscillators  of  around  25V  when  allowance  is  made  for  back- 
bias  effects. 


5. 


4.3. 


S3P  results 


The  small  signal  admittance  of  the  I'lPATT  diode  was  evaluted  over  the 
frequency  range  30  to  50GHz. 

Figure  5a  shows  an  admittance  plot  for  the  unilluminated  diode.  The 
conductance  is  seen  to  be  negative  over  the  full  frequency  range. 

Figure  5b  shows  an  admittance  plot  for  an  optically  generated  hole 

+ 

current  (corresponding  to  illumination  from  the  N  direction)  of 

-2 

50A  cm  (equivalent  to  a  current  of  1.43mA  in  a  device  of  the  area 
used  experimentally).  The  magnitude  of  the  diode  conductance  is  seen 

to  be  reduced  and  the  susceptance  changed.  Figure  5c  shows  an 

-2 

admittance  plot  for  an  optically  generated  electron  current  of  50A  cm 

+ 

(corresponding  to  illumination  from  the  P  direction).  The  change  in 
admittance  is  seen  to  be  much  greater  than  for  hole  injection, 
confirming  the  effect  of  the  higher  avalanche  ionisation  coefficient  for 
electrons  in  silicon. 

4.4.  F^P  Results 

Reference  to  equation  2  shows  that  the  oscillator  frequency  change  for 

a  given  optically  generated  current  is  expected  to  be  inversely 

proportional  to  oscillator  Q.  An  object  of  the  large  signal  modelling 

work  was  to  investigate  the  effect  of  changing  the  oscillator  Q  on 

optical  control  performance.  For  this  purpose  two  different  oscillator 

circuits  were  used.  The  first  circuit,  shown  in  Figure  6a,  is  similar  to 

that  used  in  the  previous  W-Band  work  (see  Figure  5.7  of  [1]).  It 

consists  of  an  inductor,  L  ,  selected  to  cancel  the  capacitive 

P 

susceptance  of  the  diode  at  the  oscillation  frequency  with  a  quarter- 
wave  transformer  to  match  the  diode  admittance  to  that  of  the  load. 

Close  to  resonance  the  quarter-wave  transformer  and  load  can  be 

represented  as  a  series  L-C-R  circuit,  following  the  derivation  of 
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Appendix  r where 


with  f  being  the  resonant  frequency,  Z  the  characteristic  impedance  of 
o  t 

the  transformer  and  Z^  the  load  impedance. 


Thus  the  oscillator  Q  can  be  written 


w (L  +L  ) 
P  e 


where  w  is  the  oscillation  frequency.  For  u)  =  2  n  f  L  is  series 

o  e 

resonant  with  C  and  an  equivalent  admittance  is  presented  to  the  diode 
e 

by  the  circuit  of  Figure  6b.  The  Q  for  this  circuit  is  given  by 


V  Rp 


and  is  therefore  les3  than  that  for  the  circuit  of  Figure  6a,  although  the 
admittance  presented  to  the  diode  is  the  same,  giving  equivalent 
oscillation  conditions. 
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For  oscillation  at  Ka-Band  the  following  values  were  chosen: 

R  =  1.0 

P 

L  =  0.085  nH 

P 

C  =  15pF 


Z  15 

t 

giving  R  =  4.5  and  L  =  0.049nH  at  the  design  frequency,  f  = 
e  e  o 

38GHz.  Thus  the  Q  for  the  circuit  of  Figure  6a,  from  equation  8,  is 

5.82  while  for  that  of  Figure  6b,  from  equation  9,  it  is  3-69. 

-5  2 

A  device  area  of  2.8  x  10  cm  was  used  throughout  the  modelling  work, 
corresponding  to  the  mesa  diameter  of  about  60pm  used  experimentally. 
For  the  high  Q  circuit  with  no  optically  generated  current  the  oscillation 
stabilised  at  a  frequency  of  38.708GHz  with  15-6%  voltage  modulation, 
while  for  the  low  Q  circuit  the  frequency  of  oscillation  was  39.197GHz 
with  16.2$  voltage  modulation.  The  values  are  not  identical  because  the 
quarter-wave  transformer  is  being  operated  slightly  off  its  resonant 
frequency  of  38GHz. 

Figure  7  shows  the  change  in  oscillation  frequency  for  the  two  circuits 

with  an  optically  injected  current  comprising  equal  numbers  of  electrons 

and  holes.  The  tuning  slope  for  the  low  Q  circuit  is  11.9MHz  mA 

-1 

while  that  for  the  high  Q  circuit  is  7.9MHz  mA  ,  a  ratio  of  1.51*  The 
ratio  of  the  circuit  Q  factors  is  1.58  thus  confirming  very  satisfactorily 
the  prediction  of  equation  2. 

Attempts  were  made  to  model  high  Q  circuits  representing  those  used  in 
the  experimental  work.  However,  it  was  found  that  the  circuit 


8. 


configurations  available  in  the  model  gave  rise  to  sub-harmonic 
instabilities  when  the  necessary  additional  elements  were  introduced  and 
modifications  to  the  program  would  be  necessary  to  give  a  more  general 
oscillator  circuit  modeL  Since  these  modifications  would  increase 
program  execution  time  considerably  they  have  not  been  implemented  as 
part  of  this  study. 

The  FSP  was  also  used  to  carry  out  a  more  detailed  investigation  of 

optical  tuning  transients.  Figure  8  shows  the  tuning  transient  resulting 

-2 

from  a  step  ftinction  optical  input  producing  5A  cm  in  the  low  Q 
oscillator  circuit.  As  expected  electron  injection  produces  a  larger 
tuning  effect  than  hole  injection,  owing  to  the  larger  avalanche 
Ionisation  coefficient  for  electrons  than  for  holes  in  silicon.  The  exact 
shape  of  the  initial  transient  depends  on  the  paint  in  the  oscillation 
cycle  when  the  optical  signal  is  applied.  Clearly,  tuning  is  very  fast; 
90%  of  the  frequency  change  takes  place  in  less  than  0.5ns. 

Figure  9  shows  optical  tuning  transients  for  step  function  inputs  to  the 

high  Q  oscillator  circuit.  Again  electron  injection  produces  a  much 

larger  tuning  effect.  The  curves  for  hole  injection  currents  of  5  and 
-2 

10A  cm  verify  that  the  change  in  frequency  is  proportional  to  the 
injected  current.  The  time  taken  for  90%  of  the  frequency  change  to 
take  place  is  about  0.2ns,  and  this  seems  independent  of  the  magnitude 
of  the  change. 

The  reason  for  the  difference  in  tuning  speed  between  the  two  circuits 
is  not  clear,  since  the  stored  energy  in  the  low  Q  circuit  would  be 
expected  to  be  smaller.  It  is  posable  that  the  time  at  which  the 
control  signal  has  been  applied  is  less  favourable  in  the  first  case. 
Further  examination  of  the  effect  of  applying  the  tuning  signal  at 
different  points  in  the  oscillation  cycle  is  needed  to  clarify  the 
situation.  What  can  be  asserted  with  confidence  is  that  the  tuning  is 
fast  relative  to  conventional  oscillator  control  techniques. 
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Experimental  Work 

5.1.  Device  Fabrication 

A  shift  of  emphasis  to  Ka-Band  took  place  as  part  of  the  present 
programme.  This  involved  the  fabrication  of  new  IMPATT  structures 
including  a  special  device  design  allowing  illumination  from  beneath  the 
heatsink. 

Figure  10  shows  the  basic  mesa  structure.  The  epitaxial  layer  was 

+  19  -3 

grown  by  VPE  on  an  N  substrate  of  doping  density  4.0  x  10  cm  .  An 

16  -3  + 

N  layer  doping  density  of  3-0  x  10  cm  was  used.  The  P  layer  was 

20 

produced  by  boron  diffusion  with  a  surface  concentration  of  1.0  x  10 

cm  .  A  Cr-Au  contact  was  evaporated  onto  the  P+  surface  and  the 

heatsink  plated  over  this  layer  before  mesa  etching  and  cleaving.  Mesa 

diameters  used  ranged  fV-om  50  to  65  ixm  giving  device  areas  of  from 
—5  —5  2 

1.96  x  10  to  3.3  x  10  cm  . 

+ 

In  order  to  permit  illumination  fhom  the  P  region  special  IMPATT 

structures  were  fabricated  using  a  modified  integral  heat  sink  process  to 

produce  a  25  p  m  diameter  hole  beneath  the  device.  Figure  11  is  a 

drawing  of  the  modified  device  while  Figures  12  to  15  are  scanning 

electron  micrographs  of  a  typical  experimental  device.  The  doping 

profiles  for  these  devices  were  similar  to  the  conventional  structures, 

„  -5  2 

with  a  typical  Junction  area  of  2.85  x  10  cm  . 

5.2.  Optical  Tuning  Results 

The  difficulties  caused  by  reflection  in  the  quartz  ring  package  were 
referred  to  in  [1],  For  the  present  work  250pm  quartz  stand-off  blocks 
were  used  to  permit  direct  optical  access  to  the  active  device. 

As  in  the  earlier  W-Band  work  resonant  cap  oscillator  circuits  were  used 
[8].  For  the  conventional  devices  a  50/125  pm  fibre  was  aligned 


perperi.ii  ilv  io  the  plane  of  the  Junction  using  an  XYZ  stage.  With 
this  configuration  optically  generated  currents  of  up  to  31  u  A  were 
achieved  using  a  semiconductor  laser  source  of  emission  wavelength 
850nm  and  output  power  2.9  mW.  Assuming  an  internal  quantum 
efficiency  of  unity,  which  is  a  reasonable  approximation  at  this 
wavelength,  the  coupling  efficiency  is  given  by 


where  I  is  the  optically  generated  current,  P  the  optical  power, 
so  op 

X  the  optical  wavelength,  h  Planck's  constant,  e  the  electronic  charge 

and  c  the  velocity  of  light  in  vacuo.  For  the  values  quoted  above  nc  is 

1.6?.  This  value  is  about  50?  better  than  that  obtained  using  the 

quartz  ring  package  and  W-Band  diodes  [1]  but  is  still  very  low. 


Light  emerging  from  the  cleaved  end  of  the  fibre  diverges  with  an  angle 
o 

of  about  35  .  The  structure  of  the  diode  mount  requires  the  end  of  the 

fibre  go  be  at  least  125pm  from  the  diode  mesa,  giving  a  beam  area  at 

2 

the  mesa  edge  of  12,300  urn  .  Since  the  diode  is  about  60  u  m  in 
diameter  and  5pm  high,  the  coupling  efficiency  cannot  exceed  2.4?. 
Account  also  needs  to  be  taken  of  the  reflection  loss  at  the  air/silieon 
interface.  For  normal  incidence  this  is  given  by  the  Fresnel  equation 


where  n  is  the  refractive  index  of  silicon. 

3 

Using  n  s  3.44  gives  n  =70?.  Thus  the  overall  coupling  efficiency 
s  r 

would  be  1.7?,  a  value  in  reasonable  agreement  with  the  meaatred 
results.  Any  slight  tilting  of  the  diode  mount  will  shade  part  of  the 
diode  lowering  the  coupling  efficiency  futher.  Note  that  this 
calculation  ignores  recombination  in  the  substrate  of  the  diode:  since 
the  diffusion  length  in  heavily  doped  N  type  silicon  is  of  the  order  of 


20  p  m  this  Is  a  good  approximation.  It  is  also  assumed  that  all  of  the 
light  Incident  on  the  diode  is  absorbed;  thi3  is  justified  since  the 
absorption  length  in  silicon  for  a  wavelength  of  850nm  is  about  10  urn,  a 
value  small  relative  to  the  mesa  diameter. 

Figure  16  shows  an  optical  tuning  curve  for  a  Ka-Band  oscillator  under 
edge  illumination.  The  maximum  frequency  tuning  obtained  was  40MHz, 
some  four  times  greater  than  in  the  earlier  W-Band  work  [1];  a  tenfold 
increase  in  percentage  tuning. 

The  oscillator  Q  was  measured  by  injection  Locking  to  be  about  30, 

which  with  an  estimated  voltage  modulation  of  15$  and  the  avalanche 

zone  length  estimated  in  section  4.2  gives  a  predicted  tuning  slope  of 

1.2MHz  pA  ;  this  is  in  reasonable  agreement  with  the  measured  value  of 
-1 

1.5MHz  pA  . 

Some  preliminary  measurements  on  the  modified  IMPATT  structure  were 

made.  Light  was  incident  on  the  underside  of  the  diode  through  the 

mount  shown  in  Figure  17.  Alignment  of  the  core  of  the  fibre  was  very 

difficult  to  achieve  and  the  maximum  optically  generated  current 

obtained  was  1.4  pA.  This  low  value  results  from  the  fibre  end  being 

some  100p.m  from  the  device.  Allowing  for  the  beam  divergence  from 

the  end  of  the  fibre  would  give  a  beam  area  at  the  device  of  about 
2 

9,000  p  m  :  because  of  the  limited  diameter  of  the  hole  in  the  diode 

mount  considerable  scattering  from  the  walls  of  the  hole  takes  place. 

The  hole  diameter  In  the  IMPATT  integral  heatsink  is  25  pm,  giving  an 

2 

area  available  for  illumination  of  490  pm  .  Thus  the  coupling  efficiency 
cannot  exceed  5$.  Fhom  equation  10  the  measured  coupling  efficiency 
was  about  0.1$.  Allowing  for  reflection  loss  at  the  air/sQicon  interface 
only  reduces  the  discrepancy  slightly.  It  is  thought  that  a  considerable 
improvement  in  coupling  efficiency  could  be  obtained  using  a  micro- 
positioning  system  for  fibre  alignment. 

Although  the  optically  generated  current  was  small,  significant  tuning 


12. 


was  observed.  Figure  18  shows  a  measured  tuning  characteristic.  The 

-1 

measured  tuning  slope  was  3-  1MHz  p  A  .  This  is  about  double  the  tuning 

slope  for  edge  illumination  of  similar  devices.  When  account  is  taken  of 

the  lower  bias  current,  and  hence  reduced  output  power  and  avalanche 

gain  coefficient,  F  f  X  I,  the  increased  tuning  slope  is  more  significant 
a 

stilL  The  explanation  for  this  is  probably  the  increase  in  the  proportion 

of  optically  generated  electron  current  entering  the  avalanche  zone, 

+ 

resulting  from  illumination  from  the  P  direction.  This  could  be 

investigated  by  using  shorter  wavelength  light  to  concentrate  absorption 
+ 

in  the  P  region. 


5-3-  Tuning  Speed  Measurements 

Measurements  of  optical  tuning  speed  were  carried  out  for  both  W-Band 
and  Ka-Band  oscillators.  The  laser  was  biased  above  threshold  and 
driven  with  pulses  of  Ip  s  duration  and  5ns  rise  time.  The  .change  in 
IMPATT  oscillator  frequency  was  measured  using  the  delay  line 
discriminator  system  shown  in  Figure  19. 

Figure  20  shows  the  optical  tuning  transient  for  a  W-Band  oscillator. 
The  frequency  change  of  1.4MHz  took  place  in  under  20ns.  Figure  21 
shows  the  tuning  transient  for  a  Ka  band  oscillator.  The  frequency 
change  of  5.5MHz  took  place  in  about  8ns. 

In  both  of  these  results  the  rise  time  of  the  bias  pulse  was  a  significant 
limiting  footer.  When  biased  above  threshold  the  laser  output  rise  time 
would  be  expected  to  be  under  Ins.  The  combined  discriminator  delay 
and  oscilloscope  rise  time  was  about  4.5ns  for  the  W-Band  system  and 
6ns  far  the  Ka-Band  system.  Thus  the  Ka-Band  result  is  at  the 
measuring  limit  of  the  test  set-up,  while  the  W-Band  result,  which  was 
for  a  high  Q  oscillator  circuit  (0=  200),  is  somewhat  slower.  Both 
results  are  valuable  In  that  they  confirm  that  slow  processes,  such  as 
trapping,  are  not  significant  in  the  optical  tuning  results  presented. 


13. 


5.4. 


Etching  of  Optical  Fibre 


As  discussed  in  Section  5-2  there  were  considerable  difficulties  in 
locating  the  optical  fibre  close  to  the  underside  of  the  IMPATT  structure 
so  as  to  achieve  base  illumination.  Since  the  cladding  of  the  fibre  was 
of  diameter  125  pm,  with  a  core  diameter  of  50  Pm,  considerable  scope 
exists  for  reducing  the  cladding  diameter  without  impairing  the  optical 
confinement.  The  fibre  could  then  be  located  immediately  under  the 
device  integral  heat  sin*. 

The  method  used  to  reduce  the  fibre  diameter  was  wet  chemical  etching 
with  a  1:9  mixture  of  hydrochloric  and  hydrofluoric  acid.  An  etch  rate 
of  25  pm  diaraeter/15minutes  was  observed,  with  excellent  uniformity. 
Figure  22  shows  the  transition  between  the  etched  and  unetched  fibre. 
This  technique  will  be  useful  in  further  work  on  base  illumination. 

6.  Conclusions 

A  detailed  analytic  theory  for  the  optical  control  of  IMPATT  diodes  has  been 
developed.  This  shows  that  avalanche  multiplication  of  the  optically  generated 
carriers  constitutes  an  important  gain  mechanism,  making  IMPATT  devices 
especially  suitable  for  optical  control.  For  a  typical  Ka-Band  oscillator  of  Q 
factor  30  a  1J  change  in  frequency  is  predicted  for  an  optically  generated 
current  of  leas  than  400p  A:  this  could  be  produced  by  a  ImW  laser  given  good 
coupling  efficiency. 

Calculations  of  optical  injection  locking  range  predict  that  for  a  similar  oscillator 
a  range  of  5%  should  be  achievable  with  a  peak  injected  locking  current  of  lOOpA. 


The  analytic  theory  suffers  two  main  limitations.  First,  the  power  law  avalanche 
ionisation  coefficient  model  does  not  account  for  saturation  at  high  electric 
fields;  this  makes  the  predictions  somewhat  optimistic  at  high  voltage  swings. 
Second,  the  theory  assumes  equal  avalanche  ionisation  coefficients  for  holes  and 


electr >•!'■.  T'.tj  '.t  cannot  account  for  variations  in  the  composition  of  the 
optically  generated  current  entering  the  avalanche  zone,  resulting  from  different 
illumination  configurations. 

To  remove  these  limitations  and  to  permit  studies  of  transient  effects  a 
comprehensive  computer  modal  has  been  devloped.  Among  the  important  results 
obtained  is  that  the  optical  tuning  slope  for  electron  injection  to  the  avalanche 
zone  of  a  silicon  IMPATT  is  about  four  times  that  for  hcde  injection;  a 
consequence  of  the  different  ionisation  coefficients  for  electrons  and  holes  in 
silicon.  Studies  of  the  small  signal  admittance  over  a  wide  frequency  range  have 
confirmed  the  greater  effect  of  electron  injection.  A  comparative  study  of  two 
oscillator  circuits  of  different  Q  factor  has  been  made.  This  confirmed  the  Q 
dependence  of  the  tuning  slope  predicted  by  the  analytic  theory. 

Further  transient  tuning  studies  have  been  carried  out.  These  confirm  that  sub- 
nano-second  tuning  speeds  are  readily  possible.  However,  the  reasons  for  the 
different  tuning  speed  behaviour  of  the  two  oscillator  circuits  examined  is  not 
well  understood.  Further  tests,  coupled  with  an  extended  measurement  programme 
would  be  needed  to  clarify  this  issue. 

A  major  advance  in  the  experimental  work  has  been  the  fabrication  of  a  special 

+ 

IMPATT  structure  permitting  Illumination  of  the  P  region  of  the  device. 
Preliminary  results  indicate  an  improved  tuning  slope  relative  to  edge  illumination. 
However,  the  tuning  range  was  severely  limited  owing  to  difficulties  In  bringing 
the  fibre  end  dose  enough  to  the  device  to  give  efficient  coupling.  A  wet 
etching  technique  has  been  developed  for  reducing  the  diameter  of  the  fibre  end, 
and  this,  combined  with  refinements  to  the  alignment  procedure,  would  be 
expected  to  give  much  Improved  coupling  efficiency.  An  important  area  for 
further  work  Is  the  assessment  of  the  thermal  performance  of  the  special  IMPATT 
structure. 

Removal  of  part  of  the  Integral  heat-sink  would  be  expected  to  lead  to  a 
deterioration  In  output  power  capability.  It  was  possible  to  obtain  current 
densities  sufficient  for  oscillation  in  the  experimental  devices.  However,  a 


detailed  assessment  of  power  output,  temperature  and  reliability  would  be  very 
useful. 

Work  on  edge  illumination  has  continued,  including  the  development  of  a  quartz 
stand-off  package  which  permits  direct  illumination  of  the  device.  This  has 
enabled  a  tuning  range  of  40M  Hz  to  be  obtained  at  Ka-Band,  an  improvement  of 
an  order  of  magnitude  in  percentage  tuning  over  earlier  W-Band  results  using  a 
quartz  ring  package.  There  remains  much  scope  for  further  improvement  in 
coupling  efficiency,  which  could  enhance  the  tuning  ranges  obtained  by  a  factor 
of  over  51.  Suggested  techniques  would  include  the  use  of  fibre  microlenses  to 
give  shaped  control  beams. 

Measurements  of  optical  tuning  speed  on  W  and  Ka-Band  oscillators  confirmed 
tuning  speeds  in  the  nanosecond  region.  It  was  difficult  to  correlate  the 
observed  results  with  known  oscillator  circuit  parameters,  however  the  Ka-Band 
result  of  about  8ns  was  limited  by  the  measurement  system.  Further  experiments 
with  oscillators  of  differing  Q,  operating  at  the  same  frequency,  would  be 
helpful. 

In  conclusion,  optical  tuning  of  IMPATT  oscillators  has  been  shown  to  be  a 
practical  technique  at  millimetre- wave  frequencies.  The  first  millimetre- wave 
IMPATT  structure  specially  designed  for  optical  control  has  been  fabricated  and 
demonstrated  successfully.  There  remains  a  great  deal  of  scope  for  further 
technology  development  based  on  the  work  presented  in  this  report.  Optically 
controlled  IMPATT  oscillators  having  the  tuning  performance  achieved  here  would 
be  useful  in  short  range  millimetre- wave  communication  links.  Improvements  in 
optical  coupling  efficiency  should  permit  their  use  in  phased -array  radar  and 
other  large  microwave  systems,  where  the  benefits  of  optical  control  signal 
distribution  would  be  especially  marked. 
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Figure  Captions 


Figure  la 


Variation  In  avalanche  gain  coefficient  F  (xl  with  voltage 

a  3 

swing  parameter  X. 


Figure  1b  Variation  in  avalanche  gain  coefficient  F  { Xl  with  voltage 

b  *•  3 

swing  parameter  X. 


Figure  2  Modelled  doping  profile  of  Ka-Band  IMPATT  structure. 


Figure  3  Measured  doping  profile  of  Ka-Band  IMPATT  structure. 

Figure  4  Modelled  electric  field  and  electron  current  density  profiles  for 

Ka-Band  IMPATT  diode. 


Figure  5a  Modelled  small  signal  admittance  x>f  Ka-Band  IMPATT  structure: 

urrilluminated  diode. 


Figure  5b  Modelled  small  signal  admittance  of  Ka-Band  IMPATT  structure: 

-2 

50  A  cm  hole  current. 


Figure  5c  Modelled  small  signal  admittance  of  Ka-Band  IMPATT  structure: 

-2 

50  A  cm  electron  current. 

Figure  6a  High  Q  circuit  model  for  IMPATT  oscillator  optical  tuning 

studies. 


Figure  6b  Low  Q  circuit  model  for  IMPATT  oscillator  optical  tuning  studies. 

Figure  7  Optical  tuning  characteristics  of  low  Q  and  high  Q  oacfllator 

circuits. 


Figure  8 


Optical  tuning  transients  for  low  Q  oscillator  circuit. 


Figure  9 


Optical  tuning  transients  for  high  Q  oscillator  circuit. 
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10  Cross  section  of  Ka-Band  IMPATT  structure. 

11  Drawing  of  Ka-Band  IMPATT  structure. 

12  Top  view  of  modified  Ka-Band  IMPATT  structure. 

13  Side  view  of  modified  Ka-Band  IMPATT  structure. 

14  Bottom  view  of  modified  Ka-Band  IMPATT  structure  showing  hole 

+ 

for  illumination  of  P  region. 

15  Detail  of  heat-sink  and  hole. 

16  Optical  tuning  characteristic  of  Ka-Band  IMPATT  oscillator  with 

edge  illumination  (I  =  81mA,  f  =  38.4GHz). 
do  o 

17  Mount  for  base  illumination  of  modified  IMPATT  structures. 

18  Optical  tuning  characteristic  of  modified  Ka-Band  IMPATT 

structure  with  base  illumination  (I  =  49mA,  f  =  42GHz). 

dc  o 

19  Delay  line  discriminator  system  used  for  tuning  speed 
measurements. 

20  Measurement  of  W-Band  oscillator  tuning  speed.  Upper  trace: 

laser  bias  pulse.  Lower  trace:  discriminator  output.  Horizontal 
-1 

scale  20ns  div  . 

Measurement  of  Ka-Band  oscillator  tuning  speed.  Upper  trace: 

laser  Was  pulse.  Lower  trace:  discriminator  output.  Horizontal 
-1 

scale  5ns  div  . 


Figure  21 
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Figure  22 


Etching  of  optical  fibre.  Transition  from  125pm  diameter  fibre 
to  70pm  diameter  etched  fibre.  ;(x  60  magnification). 
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FIGURE  2 


MODELLING  DOPING  PROFILE  OF  Ka-BrtND  IMPATT  STRUCTURE 
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FIGURE  A.  MODELLED  ELECTRIC  FIELD  AND  ELECTRON  CURRENT  DENSITY 
PROFILES  FOR  Ka-BAND  IMPATT  DIODE. 
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FIGURE  5a.  MODELLING  SMALL  SIGNAL  ADMITTANCE  OF  Ka-BAND  IMPATT 
STRUCTURE: 

UNILLUMINATED  DIODE 
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FIGURE  5b. 


MODELLING  SMALL  SIGNAL  ADMITTANCE  OF  Ka-BAND  IMPATT 
STRUCTURE  :  50Acm~'  HOLE  CURRENT. 
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FIGURE  5c.  MODELLING  SMALL  SIGNAL  ADMITTANCE  OF  Ka-BAND  IMPATT  STRUCTURE: 
50Acro_J  ELECTRON  CURRENT. 
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FIGURE  7. 


OPTICAL  TUNING  CHARACTERISTICS  OF  LOW  0  AND  HIGH  0 
OSCILLATOR  CIRCUITS. 
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FIGURE  g.  OPTICAL  TUNING  TRANSIENTS  FOR  HIGH  Q  OSCILLATOR  CIRCUIT 
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FIGURE  13. 


SIDE  VIEW  OF  MODIFIED  Ka-BAND  IMPATT  STRUCTURE 


FIGURE  14. 


BOTTOM  VIEW  OF  MODIFIED  Ka-BAND  IMPATT  STRUCTURE 
SHOWING  HOLE  FOR  ILLUMINATION  OF  P+  REGION. 
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FIGURE  16.  OPTICAL  TUNING  CHARACTERISTIC  OF  Ka-BAND  IMPATT 
OSCILLATOR  WITH  EDGE  ILLUMINATION. 

(  Idc  =  81mA,  f0  r  38 . AGHz ) 
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FIGURE  ?1.  MEASUREMENT  OF  K; 5-BAND  OSCILLATOR  TUNING  SPEED. 

UPPER  TRACE:  LASER  BIAS  PULSE.  LOWER  TRACE: 
DISCRIMINATOR  OUTPUT. 

HORIZONTAL  SCALE  5ns  DIV1  . 


IGURE  22. 


ETCHING  OF  OPTICAL  FIBRE.  TRANSITION  FROM  125um 
DIAMETER  FIBRE  TO  TOpim  DIAMETER  ETCHED  FIBRE. 

(X  60  MAGNIFICATION). 


A?D£NDIX  A 


General  Analytic  Theory  of  Optical  Control 


The  theory  of  optical  tuning  presented  previously  [1]  used  truncated  Bessel 
function  series  to  estimate  the  optical  tuning  capabilities  of  IMPATT  oscillators. 
As  the  voltage  swing  across  the  diode  increases  higher  order  Bessel  function 
coefficients  become  significant  and  the  model  becomes  less  accurate.  In  this 
section  more  general  expressions  are  developed  which  allow  greater  accuracy  to 
be  obtained  at  higher  voltage  swings.  Also,  opportunity  has  been  taken  to 
develop  the  theory  for  optical  injection  locking  of  IMPATT  oscillators. 

In  optical  injection  locking  the  optical  control  signal  is  intensity  modulated  at 
a  frequency  close  to  the  fhee-running  frequency  of  the  oscillator.  This  enables 
the  oscillator  frequency  to  be  locked  to  the  modulation  frequency  in  a  manner 
analogous  to  conventional  electrical  injection  locking  [2].  Whilst  optical  injection 
locking  has  yet  to  be  demonstrated  at  millimetre-wave  frequencies  because  of  the 
technological  difficulty  in  intensity  modulating  light  at  sufficiently  high 
frequencies,  recent  developments  in  optical  heterodyne  technology  [3]  suggest 
that  the  technique  may  become  of  considerable  systems  interest  in  the  future. 


2.1  IMPATT  Diode  Model 


A  general  analysis  of  the  optically  controlled  IMPATT  oscillator  is  made 
difficult  by  the  extreme  non-linearity  of  the  avalanche  multiplication  process. 
However,  if  it  is  assumed  that  all  carrier  motion  is  at  saturated  velocity,  v3, 
and  that  electrons  and  holes  have  equal  ionisation  coefficients,  a,  a  simple 
expression  for  the  avalanche  current  density,  J,  can  be  written  [H]: 


where  La  is  the  avalanche  zone  length,  xa  (-  La/vs)  is  the  avalanche  zone 
transit  time  and  Jg  is  the  total  reverse  saturation  current  density  including 
both  optically  and  thermally  generated  components: 

Js  *  J30  +  J3l  sin(w  t  +  $>)  (2) 


Here  JSQ  is  the  constant  thermally  and  optically  generated  saturation 
current  and  J3l  is  the  peak  value  of  an  optically  generated  locking  current  of 
frequency  w  and  arbitrary  phase  constant  <j>. 


To  solve  equation  ( 1 )  an  expression  for  the  electric  field  dependence  of  the 
ionisation  coefficient,  a,  is  required.  A  useful  model  is 


where  E  is  the  electric  field,  m  is  an  ionisation  non-linearity  constant  (m  *  6 
for  silicon)  and  a0  and  EQ  are  constants.  Assuming  that  the  avalanche  zone  is 
uniformly  doped,  that  the  field  in  the  drift  zone  is  much  less  than  the  peak 
field  within  the  diode  and  that  the  voltage  swing  across  the  diode  is  moderate, 
then,  following  Carroll  [5],  the  Ionisation  integral  in  equation  (1)  can  be 


evaluated  to  give 
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where  harmonic  terms  are  neglected  and  the  peak  electric  field,  Ep,  is  repre¬ 
sented  by 

Ep  -  Epo  -  Eb  ♦  Et  sin(w  t  +  6)  (5) 

for  a  diode  terminal  voltage 


v  -  VpQ  +  L(E1  ain(u  t)  -  Efa) 


(6) 


Here  Ep  is  a  back  bias  field  resulting  from  the  combined  effects  of  ionisation 
non-linearity,  reverse  saturation  current  and  carrier  space  charge,  E1  is  the 
peak  value  of  the  modulated  component  of  electric  field,  L  is  the  total  deple¬ 
tion  width  and  Vpo  the  breakdown  voltage.  Equation  (6)  assumes  the  diode  to  be 
embedded  in  an  oscillator  circuit  of  moderately  high  Q  factor,  so  that  the 
voltage  variation  will  be  nearly  sinusoidal,  and  that  the  diode  remains 
punched- through  over  the  whole  r.f.  cycle.  It  is  also  assumed  that  the  oscil¬ 
lator  is  locked  to  the  incident  modulated  optical  signal,  where  present,  to 
give  an  oscillation  frequency  of  w.  The  Injection  locking  range  can  then  be 
found  from  the  limiting  conditions  on  the  phase  angle,  $.  The  phase  angle,  6, 
is  used  to  represent  back  bias  and  carrier  space  charge  effects:  details  of  its 
evaluation  are  given  by  Cullen  and  Forrest  [6]. 


Equation  ( 4 )  is  a  linear,  first  order  equation,  and  can  be  solved  using  an 
integrating  factor  expanded  in  terms  of  modified  Bessel  functions.  Following 
integration  across  the  drift  zone  the  fundamental  component  of  the  device 


external  conduction  current  can  be  written 
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where  td  is  the  drift  zone  transit  time  and  In(X)  are  modified  Bessel  function 
coefficients  of  order  n  and  argument 
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In  deriving  equation  (7)  it  has  been  assumed  that  Jgl  <  JgQ  €  Jdc  so  that  the 
effect  of  the  optically  Injected  current  on  the  diode  bias  voltage  is  small. 


Equation  (7)  contains  three  clearly  identifiable  terms;  the  term  in 
cosCcn  t  ♦  6)  represents  the  fundamental  component  of  the  main  avalanche  cur¬ 
rent,  that  in  sin(w  t  +  6)  a  quadrature  term  that  alters  the  phasing  of  the 
avalanche  current  in  response  to  reverse  saturation  current  J30,  while  the  term 
in  costiu  t  +  <j>)  represents  the  optical  injection  locking  signal.  The  first  two 
terms  can  be  combined  by  introducing  a  modified  phase  angle,  6S: 
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where  the  avalanche  gain  factor  for  optical  tuning,  Fa{X),  is  defined  by 
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Equation  (7)  can  then  be  rewritten 
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where  Fb{X)  is  an  avalanche  gain  factor  for  optical  injection  locking: 
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Equations  (9)  to  (12)  represent  a  generalisation  of  the  earlier  analysis  [1] 
and  enable  predictions  of  optical  tuning  and  injection  locking  performance  to 
be  made  when  coupled  with  a  suitable  oscillator  circuit  model. 


The  avalanche  gain  factors  defined  by  equations  (10)  and  (12)  are  strong 
functions  of  the  peak  r.f.  voltage  swing  across  the  IMPATT  diode  through  the 
Bessel  function  argument  dependence  on  peak  modulated  electric  field,  equation 
(8).  In  evaluating  these  gain  factors  the  number  of  terms  in  the  Bessel  func¬ 
tion  series  that  need  to  be  taken  into  account  depends  on  the  argument  value, 
X.  A  typical  IMPATT  might  have 

T 3  -  0.2  Td  -  0.2  n/w0  (13) 

where  uQ  is  the  free  running  oscillator  frequency.  Thus,  for  w  *  uiQ,  m  -  6  and 
a  maximum  voltage  modulation  of  25$,  X  £  5.6.  Figure  1  shows  the  variation  in 
Fa(X}  and  F b { X }  with  X  for  n  £  1  to  4.  It  can  be  seen  that  truncating  the 
series  at  n  -  3  provides  sufficient  accuracy  for  ipost  purposes. 

For  X  -  *1.5,  corresponding  to  about  20$  voltage  modulation  Fa(X}  -  78  and 
F b l X }  -  209  showing  that  avalanche  gain  enhances  the  optical  control  sensitiv¬ 
ity  considerably. 


2.2  Oscillator  Circuit  Model 

To  obtain  predictions  of  the  optical  tuning  and  locking  ranges  a  circuit 
model  for  the  oscillator  is  required.  Figure  A.1  shows  a  simple  lumped  model  in 
which  the  IMPATT  is  represented  by  an  admittance  *i  ( -  C1  +  J  B1 )  and  the  load 
by  a  parallel  LR  circuit  of  admittance  YL  ( -  ♦  J  B^) .  The  optical  injection 
locking  signal  is  represented  by  a  current  generator,  1^,  defined  from  equation 
(11)  as 


OH) 


where  A  is  the  diode  area.  From  equation  (6)  the  r.f.  voltage  across  the  diode 


d  ( Aw) 


(21  ) 


(22) 


(23) 


This  shows  that  the  optical  tuning  slope  is  proportional  to  the  avalanche  gain 
factor,  Fa { X } ,  and  Inversely  proportional  to  the  oscillator  Q.  Comparison  of 
equation  (23)  with  the  corresponding  equation  in  [1],  equation  (4.26),  shows 
that  the  only  difference  lies  in  the  definition  of  the  factor  Fa(X).  That 
defined  by  equation  ( 4 . 27 )  of  [1]  incorporates  only  the  first  term  of  equation 
(10).  Reference  to  Figure  la  shows  that  this  introduces  a  serious  error  for  a 
voltage  swing  parameter,  X,  greater  than  4,  whereas  the  inclusion  of  an  extra 
one  or  two  terms  gives  good  accuracy  up  to  the  largest  voltage  swings  likely  to 
be  encountered  in  practice. 


2.1)  Optical  Injection  Locking  Predictions 


In  optical  injection  locking  i^  *  0.  Substituting  equations  (14),  (15), 
(17)  and  (18)  in  equation  (16),  with  the  assumption  (13)  and  for  6S  «  x/2, 
gives  the  real  part  of  equation  (16)  as 
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where  I3,  -  Jg,  A. 


Applying  the  Adler  [2]  condition,  sin  $  -  ±  1,  the  locking  range  Is  found 
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Chis  can  be  restated  in  a  form  similar  to  the  AdLer  locking  equation  by 
substituting  fhom  equation  (22)  with  equation  (17)  and  assumption  (13): 
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(26) 


Thus  the  locking  range  is  proportional  to  the  optically  injected  current  I  and 
the  avalanche  gain  factor  !  X>  ,  and  inversely  proportional  to  the  oscillator  Q. 

An  estimate  of  the  locking  range  likely  to  be  achieved  in  practice  can  be 

made  using  data  fbom  the  W-Band  IMPATT  oscillator  used  in  previous  experimental 

work  [1].  This  had  a  Q  of  200  and  an  operating  frequency  of  about  9*1  GHz  with 

a  bias  current  of  100mA.  The  voltage  swing  parameter,  X,  was  approximately 

3.5,  giving  F  1  X  1  =  40  from  Figure  1b.  Substituting  these  values  in  equation 

b 

(26)  gives  a  locking  range  dependence  on  injected  modulated  current  of  1.09GHz 
mA-1,  or  a  locking  range  of  about  22MHz  for  the  illumination  configuration  used 
in  [1].  Previous  work  on  X-Band  oscillators  [7]  gave  measured  locking  ranges  of 
a  few  megahertz,  which  would  accord  with  the  predictions  of  equation  (26)  if 
differences  in  Q  and  operating  frequency  are  taken  into  account.  It  should  be 
noted  that  equation  (26)  takes  no  account  of  locking  stability.  This  cam  be 
treated  using  the  methods  described  in  [6]. 

Figure  A.  1:  Equivalent  circuit  for  optically  controlled  IMPATT  oscillator. 
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Appendix  B-  The  Quarter  Wave  Trans  forme r. 

The  analysis  starts  from  the  equation  for  the  input  imped.ir.cr; 
of  an  arbitrarily  terminated  length  of  transmission  line 
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where  Z^.  is  the  characteristic  impedance  of  the  line,  the 
impedance  of  the  load  and  i  the  length  of  the  line.  The  propagation 
constant,  2,  i s  defined  by 

P=(V=2tI  . 12) 

c  ~T 

where  X  is  the  wavelength.  In  the  quarter  wave  transformer  t  =  X/4 
and  equation  i  reduces  to 


For  oscillator  matching  transformers  Z.  <<  Z.  .  Rationalising 
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equation  A. 3. 2.1  and  neglecting  terms  in  1.  relative  to  those  in 
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With  =  Z^/5,  and  at  a  frequency  within  20'i  of  the  resonant 
frequency  of  the  transformer,  the  error  introduced  by  equation  A. 


is  less  than  0..V.,. 


From  equation 
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where  f  is  the  frequency  at  which  e  =  >./4.  Figure  is  a 

plot  of  the  real  and  imaginary  parts  of  equation  u  as  a  function 
of  frequency. 

Real  part.  ^ 
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Where 


tt / 2  equation  u 


can  be  rewritten  as 


The  quarter  wave  transformer  and  load  can  then  be  represented  by  tne 
series  resonant,  circuit  of  Figure  B.? 
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Fig.  B.2  Series  equivalent  circuit  for  quarter  wave 

transformer  near  to  resonance. 


For  Figure  B.2 
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Close  to  series  resonance  equation  6  can  be  rewritten  as 
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Comparison  with  equation  4a  gives  the  values  of  the  equivalent 
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circuit  components  as 
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and, from  the  resonance  condition 
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Optical  control  of  W-band  impatt 
oscillators 

J  F  Singleton,  A.J.  Seeds  and  S.P.  Brunt 


hij<  i  in./  unn  i/.'./fors 

Abstract:  Results  are  presented  li>:  .in.ilvtu.jl  computer  modelling  and  experimental  •  tudics  of  optically  induced 
Irequcnvx  tuning  of  W-B.tnU  t7^  MotiH/i  impatt  oscillators  Predictions  from  an  anaklical  theory  are  com¬ 
pared  with  those  from  a  large-signal  linie-domain  computet  nunfel.  and  the  effect  of  unequal  electron  and  hole 
ionisation  coefficients  on  the  tuning  performance  is  illustrated  Experimental  measurements  are  presented 
0-  >vttng  an  optical  tuning  range  ol  about  l * »  MU/  for  an  opttcaik  generated  current  of  20  ft  A.  and  comparisons 
are  made  with  the  theoretical  predictions 


1  Introduction 

The  impatt  diode  is  currently  the  onh  solid-state  source 
capable  of  pro\iding  the  transmitter  power  requirements 
in  small,  millimetre-wave  radar  systems  In  many  applica¬ 
tions.  particularly  in  KMC  W  radars,  il  is  necessary  to  tune 
the  transmitter  frequency  electronically,  however  such 
control  of  impatt  diode  oscillators  is  found  to  be  very  diffi¬ 
cult  to  implement  because  of  the  power  impedance  levels 
involved. 

Tuning  of  impatt  oscillators  in  practical  systems  is  nor¬ 
mally  achieved  by  using  a  varactor-controlled  Gunn  oscil¬ 
lator  (VCO)  to  injection  lock  the  impatt.  This  system 
works  very  well  but  has  the  disadvantage  of  an  extra  oscil¬ 
lator  cavity  which  may  cause  problems  where  size  and 
weight  of  the  system  are  prime  considerations 

A  novel  alternative  approach  is  to  use  an  optical  signal 
to  generate  carriers  in  the  impatt  by  the  photo-electric 
effect.  These  carriers  drift  into  the  avalanche  /one  of  the 
diode  and  alter  the  phasing  of  the  avalanche  cycle,  thereby 
changing  the  oscillation  frequency. 

Optical  frequency  tuning  effects  have  previously  been 
observed  in  Trapatts  ( t  f.  bipolar  transistors  [2]  and 
MESFKTS  [3]  Early  work  on  the  optical  tuning  of 
impatts  was  carried  out  predominantly  at  X-Band  fre¬ 
quencies  [4.  5].  This  paper  reports  work  to  extend  the 
technique  to  W-Band  (75  1 10  GHz). 


2  Analytic  theory  of  optical  frequency  tuning  in 
impatt  oscillators 


The  physical  mechanism  underlying  optical  frequency 
tuning  is  the  generation  of  holc-clectron  pairs  in  the  impatt 
by  the  incident  light,  with  consequent  alteration  in  the 
timing  of  the  avalanche  cycle.  Il  is  therefore  possible  to 
calculate  the  change  in  oscillation  frequency  using  a  modi¬ 
fied  Read  [6]  approach.  Making  the  simplifying  assump¬ 
tions  of  equal  carrier  velocities  and  avalanche  ionisation 
coefficients  for  holes  and  electrons,  the  total  avalanche 
current  is  given  by  [6] 


xa  dJ 
2  dt 


+  */, 


(I) 


where  Jt  is  the  total  reverse  saturation  current  density 
including  both  thermal  and  optical  carrier  generation,  n 
the  avalanche  ionisation  coefficient.  wa  the  avalanche  zone 


width  and  r „  the  avalanche  zone  transit  time.  Integration 
of  eqn.  1  using  a  power  law  approximation  for  the  depen¬ 
dence  of  the  ionisation  coefficient  on  electric  field,  together 
with  a  transit-time  model  for  the  drift  zone,  enables  the 
change  in  oscillator  frequency  for  a  small  optically  injected 
current.  to  be  written  as 
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where  /,.  IS  the  diode  bias  current.  Q  the  oscillator  external 
(/-factor  and  F,J.v}  an  avalanche  gain  factor  defined  by 
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Here  are  modified  Bessel  function  coefficients  of 

order  n  and  argument 


(m  +  m,  , 

Xl.  ru  tp.,  " 

where  /„  is  the  oscillation  frequency.  V,  the  peak  RF 
voltage  across  the  diode.  the  peak  electric  field  at 
breakdown,  w  the  total  depletion  width  and  m  the  ionis¬ 
ation  non-linearity  coefficient.  F.{. x}  is  a  strong  function  of 
the  RK  voltage  swing  though  the  K,  dependence  of  the 
Bessel  function  argument  x.  For  simplicity,  eqn.  3  neglects 
Bessel  coefficients  of  order  greater  than  2,  a  reasonable 
approximation  for  small  to  moderate  voltage  swings 
t  v  S  4). 

The  relationship  between  the  optically  generated 
current  and  the  incident  optical  power  is  dependent  on  the 
device  structure  and  illumination  geometry.  Assuming  that 
all  of  the  light  is  absorbed  within  the  depiction  layer  the 
optically  generated  current  is  given  by 


n,  ex  Pop 
he 


(5) 


where  is  the  incident  optical  power,  ij,  the  quantum 
efficiency,  a  the  optical  wavelength,  e  the  electronic  charge. 
h  Planck's  constant  and  c  the  velocity  of  light  in  vacuo. 
For  /.  =  850  nm  and  n,  =  0.9,  f„/P„,  =  0.62  A  W" '.  In 
practice,  difficulties  in  securing  efficient  coupling  between 
the  light  source  and  the  impatt  lead  to  much  smaller 
values. 
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3  Computer  modelling  studies 

The  analytical  theory  gives  a  very  useful  estimate  of  the 
magnitude  of  optical  tuning  effects,  but  it  is  limited  in 
accuracy  as  a  result  of  the  assumption  of  equal  ionisation 

14<) 


%.iH.‘tTiueni’»  tor  electron's  and  holes  and.  in  any  ease,  cannot 
he  j'eJ  to  study  transient  ejects  tn  devices. 

\  iaree-sienal  time-domain  computer  model  of  the  opti¬ 
mally  controlled  inipatt  has  therefore  been  developed  from 
the  earlier  avalanche  diode  modelling  work  of  Blakey  et  ai 
S']  I  he  model  is  based  on  the  solution  in  one  spatial 
dimension  and  time  of  the  drift -diffusion  approximation  to 
Boltzmann's  transport  equation,  the  equations  for  electron 
and  hole  continuity,  and  (iauss*  theorem  within  the  semi¬ 
conductor  material  of  the  device  This  approach  enables 
spatially  varying  optical  carrier  generation  to  be  taken 
into  account 

Modeling  studies  were  performed  for  a  W-Band, 
silicon,  single  drift  structure  having  a  doping  profile 
similar  to  that  of  the  devices  used  in  the  experimental 
work  \  simple  oscillator  circuit,  shown  in  Fig.  1,  was 


Fig  1  CKi  ir/dfi'r  egunaUnt  ciri  uit  lor  modelling  studies 


coupled  to  the  diode  to  permit  predictions  of  optical  fre¬ 
quency  tuning.  Lp  represents  the  diode  bond  wire  induc¬ 
tance.  Rp  the  diode  parasitic  resistance.  R,  the  load 
resistance.  (  the  IX’  block  capacitance  and  Z,  a  quarter- 
wavelength  transformer  If  Lp  is  selected  to  resonate  with 
the  capacitive  susceptance  of  the  diode,  then  (he  imped¬ 
ance  it  plane  X  X  .  looking  towards  the  diode,  is  real 
and  negative  The  amplitude  of  the  oscillation  is  then 
determined  by  selecting  /,  to  match  RL  to  the  impedance 
at  \  X  for  the  required  operating  point.  With  values  of 
4  pF  for  (  p .  1  O  for  Rp.  and  a  bias  current  of  140  mA.  the 
model  gave  stable  oscillation  frequency  of  94  GFfz.  The 
modelled  efficiency  was  3.6%.  which  is  a  typical  value  for  a 
W-Band  single-drift  silicon  device 

The  computer  model  can  also  be  used  to  study  transient 
effects,  and  Fig  2  shows  the  tuning  transient  resulting 
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from  u  slop  change  in  optically  injected  electron  current 
from  t)  to  2  A  cm  :  It  is  seen  [hat  Ul)  of  the  frequency 
change  was  aehtesed  in  55  ps.  sslneli  compares  sets 
favourably  with  the  speed  of  eonsenttonal  electronic 
tuning  techniques  such  as  varactor  tuning  and  mas  allord 
novel  systems  possibilities. 

Two  different  optical  carrier  generation  configurations 
were  investigated  in  the  study.  In  the  first,  electrons  were 
injected  into  the  depletion  /one  from  the  p'  region, 
whereas  in  the  second,  holes  were  injected  from  the  n ' 
region.  These  correspond  to  illumination  from  opposite 
ends  of  the  impatt  structure.  Tig.  3  shows  the  optical 


Fig.  3  Companion  between  computer  model  and  analytic  theory  for 
optical  tuning  characteristics 
Bus  current  density  IS  3  KA  cm' 


tuning  characteristics  obtained  for  optically  injected 
current  densities  of  up  to  2  A/cm2,  corresponding  to  a 
current  of  10  pA  in  the  experimental  diode.  The  change  in 
frequency  with  injected  current  is  seen  to  be  linear,  as  pre¬ 
dicted  by  eqn.  2.  A  more  detailed  comparison  with  the 
analytic  theory  can  be  made  by  substituting  appropriate 
values  into  eqns.  2-4.  For  m  -  6,  =  0.1  1  and  Q  =  1.6, 

with  the  modelled  bias  current  and  voltage  modulation, 
A/0 //„  =  8.8  MHz  jjA  This  slope  is  plotted  as  a  solid 
line  in  Fig.  3.  The  modelled  values  for  electron  and  hole 
tuning  lie  above  and  below  this  line,  respectively.  This  is 
because  the  analytic  theory  assumes  equal  avalanche  ionis¬ 
ation  coefficients  for  electrons  and  holes,  whereas  in  silicon 
those  for  electrons  are  much  greater  than  those  for  holes. 
Thus  electron  injection  produces  a  much  greater  tuning 
effect  than  hole  injection.  With  this  restriction,  agreement 
between  the  analytic  theory  and  the  modelling  results  is 
seen  to  be  good. 

4  Experimental  work 

The  devices  used  were  single-drift  silicon  impatts  mounted 
in  quartz  ring  packages  as  shown  in  Fig.  4.  All  measure¬ 
ments  were  carried  out  using  the  well  known  resonant  cap 
circuit  design  [8]  shown  in  Fig.  5A.  At  a  current  of 
120  mA,  the  oulpul  power  obtained  was  55  mW  and  the 
efficiency  was  3.1%,  which  is  in  good  agreement  with  the 
computer  model  result  of  3.6%,  bearing  in  mind  that  the 
computer  model  neglects  circuit  losses. 

A  GaAs/GaAIAs  laser  emitting  at  850  nm  was  used 
as  the  light  source  for  the  optical  tuning  experiments. 
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The  laser  was  supplied  with  an  optical  fibre  output  of 
50  125  jim  graded  index  fibre  which  was  coupled  to  the 
edge  of  the  active  device  through  a  small  hole  in  the  oscil¬ 
lator  cavity  wall.  The  tibre  was  brought  into  very  close 
proximity  with  the  quartz  ring  of  the  diode  package,  and 
its  position  was  then  optimised  for  maximum  tuning  effect 
using  a  micromanipulator.  Fig.  5B  gives  details  of  the  con¬ 
figuration  used.  The  laser  package  also  contained  a 
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Vcriic»|  section  through  centre  of  diode 

monitor  photodiode  which  was  used,  after  calibration,  to 
determine  the  laser  output  power. 

Using  this  arrangement,  optically  generated  reverse 
saturation  currents  of  up  to  20.5  pA  were  obtained  for  an 
optical  power,  measured  at  the  fibre  end,  of  3  mW.  From 
eqn.  5  this  power  should  be  sufficient  to  generate  a  current 
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of  about  1.9  mA  for  a  quantum  efficiency  of  0.9.  The  low 
measured  coupling  efficiency  of  1.1".,  is  due  partly  to  ref¬ 
raction  in  the  quartz  ring  package  and  partly  to  the  very 
small  height  of  the  impatt  structure,  around  3  ^m,  com¬ 
pared  with  the  fibre  core  diameter  of  50  /im. 

Fig.  6.  shows  a  typical  spectrum  analyser  display  for  an 


Fig.  6  Optical  tuning  of  V.40  Mil:  at  a  centre  frequency  of  V/ .H26  GHz 
Frequency  span  20  MM/,  optically  generated  current  20.3  /iA 

oscillator  undergoing  optical  tuning.  In  this  example  the 
bias  current  was  100  mA.  giving  a  free-running  frequency 
of  91.83  GHz.  which  is  represented  by  the  peak  at  the 
centre  of  the  display.  The  right-hand  peak  is  the  output  of 
the  oscillator  with  an  optically  generated  current  of 
20.3  (/A.  The  frequency  shift  obtained  was  9.4  MHz. 

Fig.  7  is  a  plot  of  oscillation  shift  against  optically  gen¬ 
erated  current  for  another  impatt  sample  at  a  bias  current 
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Fig.  7  Frequency  tuning  against  optically  generated  current 

of  120  mA.  The  shift  is  proportional  to  current  as  predict¬ 
ed  by  eqn,  2.  In  order  to  compare  the  tuning  slope  with 
theory  it  is  necessary  to  know  the  Q-factor  of  the  oscil¬ 
lator.  This  was  measured  by  injection  locking  [9].  From 
measurements  of  the  injected  power,  the  impatt  oscillator 
output  power  and  the  locking  range  the  Q  was  determined 
from 


where  (,)„  is  the  free  running  impatt  frequency.  2Acu,  is  the 
total  locking  range.  P,  is  the  injected  power  and  P„  is  the 
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impalt  output  power.  Measurements  on  the  oscillator  of 
fig  7  gave  a  Q  of  about  200.  From  the  tuning  character¬ 
istic.  Fig.  7.  if| A/ 1  if/,,,  =  225  KHz.fjA  *.  F„|.x)  was  esti¬ 
mated  from  the  measured  output  power  to  be  27,  giving  a 
calculated  tuning  slope  of  253  KHz./<A  1  from  eqn.  2. 
There  is  thus  good  agreement  between  predicted  and  mea¬ 
sured  results,  particularly  when  bearing  in  mind  that  the 
theory  assumes  equal  ionisation  coefficients  for  holes  and 
electrons. 


5  Discussion 

Optical  tuning  has  been  demonstrated  for  a  W-Band 
impatt  oscillator.  The  maximum  frequency  shift  obtained 
was  9  4  MFfz  using  a  3  mW  semiconductor  laser.  This 
figure  is  very  small  and  can  be  explained  by  the  very  low 
coupling  efficiency  (1.1%)  of  light  from  the  laser  output  to 
the  impatt  device  due  to  the  illumination  geometry  used.  If 
this  coupling  efficiency  could  be  increased  to  65%,  which 
is  a  typical  value  for  high  speed  photodiodes,  the  tuning 
range  should  increase  to  over  600  MFlz,  which  would  be 
useful  in  many  systems  applications. 

Methods  of  improving  coupling  include  the  replacement 
of  the  quartz  ring  package  with  one  using  quartz  standoffs 
and  the  use  of  lenses  or  mirrors  to  concentrate  the  incident 
light.  A  more  promising  approach,  however,  is  to  illumi¬ 
nate  the  top  or  bottom  surface  of  the  device  rather  than 
the  edge,  as  has  already  been  demonstrated  at  lower  fre¬ 
quencies  by  Vyas  el  at.  [10]  and  Chiu  and  Freyer  [II], 
This  involves  the  fabrication  of  special  millimetre-wave 
impatt  structures  and  forms  the  basis  of  continuing  work 
on  this  project. 

Experimental  results  showed  a  linear  relationship 
between  frequency  shift  and  optically  injected  current  in 
agreement  with  the  analytical  and  computer  models  devel¬ 
oped  during  this  study.  Measured  tuning  slopes  were  also 
in  good  agreement  with  those  predicted  theoretically, 
although  the  tuning  slope  of  the  oscillator  used  in  the  com¬ 
puter  modelling  studies  was  much  larger  than  that  for  the 
experimental  oscillator,  since  a  Q  of  1.5  was  used  com¬ 
pared  with  a  Q  of  around  200  for  the  practical  oscillator. 

Tne  simple  analytical  theory  has  been  shown  to  give 
good  agreement  with  a  large-signal  time-domain  computer 
model.  The  computer  model  has  shown  that  electron  injec¬ 
tion  gives  a  considerably  greater  tuning  effect  than  hole 
injection,  which  is  an  important  consideration  when  struc¬ 
tures  are  to  be  fabricated  to  enable  the  top  or  bottom 
device  surfaces  to  be  illuminated. 


An  initial  investigation  of  the  tuning-speed  modelling 
has  shown  that  very  fast  tuning  is  possible,  which  is  ohvi- 
ouslv  a  key  factor  in  systems  applicability  It  is  expected 
that  the  oscillator  (7  would  be  the  fundamental  limitation 
in  this  regard. 

The  main  area  for  development  of  the  computer  model 
is  the  improvement  of  the  oscillator  circuit  from  the  very 
simple  one  of  Fig.  I  to  one  which  would  model  the  actual 
configuration  used.  Fig  5A. 

Optical  tuning  provides  a  technique  for  controlling  the 
frequency  of  impatt  oscillators  which  offers  electrical  iso¬ 
lation  and  ease  of  control  signal  distribution  through 
optical  fibres.  With  further  development  this  technique  will 
prove  to  be  attractive  for  a  number  of  millimetre-wave 
radar  and  communication  system  applications 
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Abstract— ThU  paper  describes  methods  of  controlling  the  frequency 
of  IMPATT  oscillators  using  optical  rather  than  electrical  signals.  An¬ 
alytic  theories  of  optical  tuning  and  injection  locking  are  presented. 
Results  from  a  comprehensive  large  signal  computer  model  of  the  op¬ 
tically  controlled  IMPATT  oscillator  are  given,  illustrating  the  Impor¬ 
tance  of  the  composition  of  the  optically  generated  current  on  the  op¬ 
tical  control  performance  and  demonstrating  the  capability  of  optical 
control  for  rapid  frequency  tuning.  Finally,  recent  experimental  work 
on  optical  tuning  effects  in  W-Band  IMPATT  oscillators  is  presented. 


I.  Introduction 

THE  EASE  of  transmission  of  high  rate  intensity  mod¬ 
ulated  optical  signals  through  single  mode  optical 
fibres  suggests  tlteir  use  for  signal  distribution  in  complex 
microwave  systems  such  as  phased-array  radars  [1],  [2], 
For  active  phased  arrays,  where  each  element  has  its  own 
microwave  power  source,  considerable  simplification  re¬ 
sults  from  direct  optical  control  of  the  active  device. 

Three  distinct  modes  of  control  can  be  identified.  In  the 
first  the  optical  signal  alters  the  negative  resistance  prop¬ 
erties  of  the  active  device  so  as  to  drive  it  into  or  out  of 
oscillation:  this  is  known  as  optical  switching.  The  sec¬ 
ond  mode  uses  an  optical  signal  of  lower  intensity  to  in¬ 
duce  optical  tuning  of  the  oscillator  frequency.  Finally, 
illumination  with  an  optical  signal,  intensity  modulated  at 
a  frequency  close  to  the  free-running  oscillator  frequency, 
one  of  its  harmonics  or  subharmonics,  can  produce  optical 
injection  locking. 

Optical  control  of  a  variety  of  microwave  devices  has 
been  demonstrated.  Yen  and  Bamoski  [3]  have  investi¬ 
gated  optical  switching  and  subharmonic  injection  lock¬ 
ing  in  bipolar  transistor  oscillators.  Optical  tuning  of  the 
oscillation  frequency  in  MESFET  oscillators  has  been 
demonstrated  by  Sun  el  at.  [4]  as  has  optical  injection 
locking  by  Salles  and  Forrest  [5].  Turning  to  two  terminal 
devices,  optical  tuning  and  switching  of  TRAPATT  os- 
cillators  operating  at  frequencies  up  to  1 .5  GHz  has  been 
reported  by  Kiehl  [6]. 

At  the  higher  microwave  frequencies  the  IMPATT 
diode  oscillator  is  one  of  the  most  powerful  solid-state 
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sources  [7]  However,  the  power  and  impedance  levels 
involved  make  tuning  by  conventional  means,  such  as  YIG 
spheres  or  varactor  diodes,  extremely  difficult.  This  has 
led  to  considerable  work  on  the  optical  control  of  IM¬ 
PATT  diodes,  in  which  optical  switching  [8],  [9],  tuning 
[10],  [II]  and  injection  locking  [12],  [13]  have  al!  been 
demonstrated.  In  addition,  optical  carrier  injection  can  as¬ 
sist  in  stabilizing  the  level  of  reverse  saturation  current  in 
the  device,  leading  to  an  improvement  in  oscillator  noise 
(141,  [15], 

The  object  of  this  paper  is  to  present  analytical  and 
computer  modelling  studies  of  optical  tuning  and  injec¬ 
tion  locking  in  IMPATT  oscillators,  together  with  recent 
experimental  measurements  of  optical  tuning  in  W-Band 
IMPATT  oscillators. 


II.  Analytic  Theory  of  Optical  Control 
A.  IMPATT  Diode  Model 

Illumination  of  an  IMPATT  diode  with  light  having  a 
photon  energy  greater  than  the  band  gap  of  the  semicon¬ 
ductor  material  from  which  it  is  fabricated  results  in  the 
generation  of  electron-hole  pairs  within  the  device.  These 
carriers  alter  the  timing  of  the  avalanche  cycle,  thus  pro¬ 
ducing  a  change  in  oscillation  frequency.  If  the  intensity 
of  the  incident  optical  signal  is  modulated  a  varying  cur¬ 
rent  is  generated  in  the  diode  which,  following  avalanche 
multiplication,  can  be  used  to  injection  lock  the  IMPATT 
oscillator. 

A  general  analysis  of  the  optically  controlled  IMPATT 
oscillator  is  made  difficult  by  the  extreme  nonlinearity  of 
the  avalanche  multiplication  process.  However,  if  it  is  as¬ 
sumed  that  all  carrier  motion  is  at  saturated  velocity  v,, 
snd  that  electrons  and  holes  have  equal  ionization  coeffi¬ 
cients  a.  a  simple  expression  for  the  avalanche  current 
density  /  can  be  written  (16): 


I .dJ 
2  dt 


(1) 


where  L,  is  the  avalanche  zone  length,  r,  (  -  L,/v,)  is 
the  avalanche  zone  transit  time  and  /,  is  the  total  reverse 
saturation  current  density  including  both  optically  and 
thermally  generated  components: 

J,  “  Ao  +  J,\  »in  (“*  +  6)-  (2) 

Here  J, 0  is  the  constant  thermally  and  optically  gener¬ 
ated  saturation  current  and  /, ,  is  the  peak  value  of  an  op¬ 
tically  generated  locking  current  of  frequency  a  and  ar¬ 
bitrary  phase  constant  d- 


0733-8724/87 /0300-0403W  1.00  ©  1987  IEEE 


404 


HR  RNAL  OF  I.IOHIRAVF  I  M  MNOl  <  X.'.  VOt  IT'  SO  I  MARCH  148' 


To  solve  (1),  an  expression  for  the  electric  field  depen¬ 
dence  of  the  ionization  coefficient  a  is  required  A  useful 
model  is 

(3) 

where  £  is  the  electric  field,  m  is  an  ionization  nonline¬ 
arity  constant  (m  «  6  for  silicon )  and  a0  and  E0  are  con¬ 
stants.  Assuming  that  the  avalanche  zone  is  uniformly 
doped,  that  the  field  in  the  drift  zone  is  much  less  than  the 
peak  field  within  the  diode  and  that  the  voltage  swing 
across  the  diode  is  moderate,  then,  following  Carroll  (17], 
the  ionization  integral  in  (1)  can  be  evaluated  to  give 

dj  =  2  (m  +  I)  Tm  / 
dt  r„  L4  \E„ 0)  E ^ 

+  -  sin  (ui /  +  6)  +  —  (4) 

£,0  • 

where  harmonic  terms  are  neglected  and  the  peak  electric 
field  Ep  is  represented  by 

Ep  =  E,#  -  Eb  +  E,  sin  (car  +  6)  (5) 

for  a  diode  terminal  voltage 

v  =  yp o  +  L(EI  sin  (uf)  -  Eb).  (6) 

Here  Eb  is  a  back  bias  field  resulting  from  the  combined 
effects  of  ionization  nonlinearity,  reverse  saturation  cur¬ 
rent  and  carrier  space  charge  X,  E ,  is  the  peak  value  of 
the  modulated  component  of  electric  field,  L  is  the  total 
depletion  width,  and  V#  the  breakdown  voltage.  Equation 
(6)  assumes  the  diode  to  be  embedded  in  an  oscillator  cir¬ 
cuit  of  moderately  high  Q  factor,  so  that  the  voltage  vari¬ 
ation  will  be  neariy  sinusoidal,  and  that  the  diode  remains 
punched-through  over  the  whole  RF  cycle.  It  is  also  as¬ 
sumed  that  the  oscillator  is  locked  to  the  incident  modu¬ 
lated  optical  signal,  where  present,  to  give  an  oscillation 
frequency  of  u.  The  injection  locking  range  can  then  be 
found  from  the  limiting  conditions  on  the  phase  angle  <h. 
The  phase  angle  i  is  used  to  represent  back  bias  and  ear¬ 
ner  space  charge  effects.'  details  of  its  evaluation  are  given 
by  Cullen  and  Forrest  [18]. 

Equation  (4)  is  a  linear,  first  order  equation,  and  can  be 
solved  using  an  integrating  factor  expanded  in  terms  of 
modified  Bessel  functions.  Following  integration  across 
the  drift  zone  the  fundamental  component  of  the  device 
external  conduction  current  can  be  written 


where  rd  is  the  drift  zone  transit  time  and  /„ )  X  )  are  mod 
ified  Bessel  function  coefficients  of  order  n  and  argument 

2  (m  +  1  )  E, 

X  =  - - (81 

wraEpo 

In  deriving  (7)  it  has  been  assumed  that  7,,  «  Jb<  « 
Jj,  so  that  the  effect  of  the  optically  injected  current  on 
the  diode  bias  voltage  is  small. 

Equation  (7)  contains  three  clearly  identifiable  terms; 
the  term  in  cos  (ut  +  S)  represents  the  fundamental  com¬ 
ponent  of  the  main  avalanche  current,  that  in  sin  (u it  + 
6)  a  quadrature  term  that  alters  the  phasing  of  the  ava¬ 
lanche  current  in  response  to  reverse  saturation  current 
Ao.  while  the  term  in  cos  (wt  +  <t>)  represents  the  optical 
injection  locking  signal.  The  first  two  terms  can  be  com¬ 
bined  by  introducing  a  modified  phase  angle  6,: 


6,  =  6  +  sin 


where  the  avalanche  gain  factor  for  optical  tuning  F„  {  X  } 
is  defined  by 


Fa{X)  =  ll{X)  + 


/,{*} 

(-l)’(2n  4  I )  /,  1 1{  Jf )  /,,{*} 
n(n  +  1) 


Equation  (7)  can  then  be  rewritten 


A  =  -2  A, 


•  “.(«(»- j) +  «.) 


Fb{X }  cos  (u(r  -  y)  +  *) 


where  Fb  {  X  }  is  cn  avalanche  gain  factor  for  optical  in¬ 
jection  locking: 


c 
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/j{X}  +  /,{*}  /,{*] 

2 


2  f 

»• 2 ("  -  1 ) ("  +  1 ) 


Equations  (9)— ( 12)  represent  a  generalization  of  an  earlier 
analysis  [19j  and  enable  predictions  of  optical  tuning  and 
injection  locking  performaixe  to  be  made  when  coupled 
with  a  suitable  oscillator  circuit  model . 

The  avalanche  gain  factors  defined  by  ( 10)  and  (12)  are 
strong  functions  of  the  peak  RF  voltage  swing  across  the 
IMPATT  diode  through  the  Bessel  function  argument  de¬ 
pendence  on  peak  modulated  electric  held,  equation  (8). 
In  evaluating  these  gain  factors  the  number  of  terms  in  the 
Bessel  function  series  that  need  to  be  taken  into  account 
depends  on  the  argument  value  X  A  typical  IMPATT 
might  have 

ra  =  0.2  Tj  -  O.2ir/u)0  (13) 

where  w0  is  the  free  running  oscillator  frequency.  Thus, 
for  co  “  oo0,  m  =  6  and  a  maximum  voltage  modulation 
of  25  percent,  X  s  5.6.  Fig  1  shows  the  variation  in 
Fa  {  X  }  and  F„  { X  }  with  X  for  n  <  I  to  4  It  can  be  seen 
that  truncating  the  senes  at  n  =  3  provides  sufficient  ac¬ 
curacy  for  most  purposes. 

ForX  =  4.5,  corresponding  to  about  20-percent  voltage 
modulation  F,  {  X  )  =78  and  Fb  {  X  }  =  209  showing  that 
avalanche  gain  enhances  the  optical  control  sensitivity 
considerably. 


B.  Oscillator  Circuit  Model 
To  obtain  predictions  of  the  optical  tuning  and  locking 
ranges  a  circuit  model  for  the  oscillator  is  required.  Fig. 
2  shows  a  simple  lumped  model  in  which  the  IMPATT  is 
represented  by  an  admittance  T,  (  =  G,  +  jB, )  and  the 
load  by  a  parallel  LR  circuit  of  admittance  Y,  (  =  GL  + 
jBt).  The  optical  injection  locking  signal  is  represented 
by  a  current  generator  i,  defined  from  equation  (11)  as 


where  A  is  the  diode  area.  From  (6)  the  RF  voltage  across 
the  diode  is  given  by 

t'l  =  ~iv i  expy(wr)  ( 15) 

where  P,  =  E,L  Thus  the  circuit  equation  for  Fig.  2  is 
i,  =  (>',  +  Y, )  t>,  (16) 


C.  Optical  Tuning  Predictions 
If  there  is  no  injected  locking  signal  GL  and  B,  are  ob¬ 
tained  from  (11),  (15).  and  (16),  setting  J,t  =  0  and  add- 
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mg  the  diode  depletion  susceptance  ujCj  to  the  active  sus- 
ceptance 


4  sin 


G ,  = 


f)i,M 


^Tj/„{x} 


/  WTj 

V  2 


5‘) 


07) 


e, 


4 ldc  sin 


Mf,) 

2 


^Tdl0{X) 


(  —d 
V  2' 


5 •) 


-  uiCd 


(  18) 


where  ldc  =  JdcA.  Dividing  (18)  by  (17)  and  defining  B, 
as  Bl  =  -  l/(uLi)  gives 

u»Q  -  ^  =  G,  cot  ^  -  8,)  ■  ( 19) 


The  tuning  sensitivity  can  be  obtained"  by  assuming  that 
ui  =  u)0  +  Aw,  where  w0  =  ( LLCd )  ~°  5.  Then  for  6,  « 
r/2  substitution  of  (9)  in  ( 19)  with  assumption  ( 1 3)  yields 


Au> 


<h 

2  Cd 


b  +  2 


r„ 


(20) 


where  =  JrfA.  The  optical  frequency  tuning  slope  is 
then 


jjAw)  =  G,F,,{X) 
<^i0  Gdldcu0ta 
Defining  oscillator  Q  factor  as 


(21)  can  be  rewritten 

<*(Aw)  =  Fa{X} 
dl*i  Ql**' 


(21) 


(22) 


(23) 


This  shows  that  the  optical  tuning  slope  is  proportional  to 
the  avalanche  gain  factor  F.  { X }  and  inversely  propor¬ 
tional  to  the  oscillator  Q 


D  Optical  Injection  Locking  Predictions 


In  optical  injection  locking,  ij  *■  0.  Substituting  (14), 
(15),  (17),  and  (18)  in  (16),  with  the  assumption  (13)  and 
for  6,  «  n /2,  gives  the  real  part  of  (16)  as 


4/„f»{X) 

*unr.Vl 


sin  $  =  mCd  - 


I 

uLl 


(24) 


where  I,,  =  J,tA . 

Applying  the  Adler  [20]  condition,  sin  <t> 
locking  range  is  found  to  be 

2  Aw,  _  *'.*1*1. 


±  1 ,  the 


(25) 


This  can  be  restated  in  a  form  similar  to  the  Adler  locking 
equation  by  substituting  from  (22)  with  (17)  and  assump- 


VOl. 


non  (13): 


2  Aw,  =  5  {  X  } 


/p  {  jf  }  Wp/,, 

/i  {X}Qld, 


(26) 


Thus  the  locking  range  is  proportional  to  the  optically  in¬ 
jected  current  /,,  and  the  avalanche  gain  factor  Fb{X  }, 
and  inversely  proportional  to  the  oscillator  Q. 


III.  Computer  Modeling  Studies 

The  analytic  theory  presented  in  the  previous  section 
requires  a  number  of  approximations  to  produce  tractable 
results.  Of  particular  importance  is  the  assumption  of 
equal  ionization  coefficients  for  electrons  and  holes.  In 
silicon  the  coefficient  for  electrons  is  an  order  of  magni¬ 
tude  higher  than  that  for  holes,  over  a  wide  rangeof  elec¬ 
tric  field.  The  optical  tuning  behavior  is  thus  considerably 
dependent  on  the  composition  of  the  optically  generated 
current  entering  the  avalanche  zone,  and  hence  on  the  il¬ 
lumination  configuration  [21]. 

To  investigate  these  effect  a  large  signal  time-domain 
computer  model  for  the  optically  controlled  1MPATT  has 
been  developed.  This  is  based  on  the  numerical  methods 
described  by  Blakey  er  al.  [22],  while  the  modifications 
required  to  model  optical  carrier  generation  effects  have 
been  detailed  elsewhere  [23],  The  model  uses  an  explicit 
solution  in  one  spatial  dimension  and  time,  with  the  active 
region  of  the  diode  divided  into,  typically,  100  mesh- 
points.  Output  from  the  model  is  in  terms  of  terminal  cur¬ 
rent  and  voltage  as  a  function  of  time. 

A.  Large  Signal  Admittance  Studies 

By  coupling  the  diode  to  an  RF  voltage  source  plots  of 
large  signal  admittance  against  frequency  and  peak  RF 
voltage  can  be  obtained.  Fig.  3  shows  typical  results  for 
an  X-Band  p*-n-n*  silicon  device,  at  a  bias  current  of 
600-A  cm  .  Carrier  generation  in  the  p+  and  n*  regions 
has  been  neglected  so  that  the  optically  generated  charge 
comprises  electron-hole  pairs  generated  in  the  n  region. 
An  exponential  optical  absorption  profile  is  assumed,  with 
an  absorption  coefficient  of  6  x  105  m"1,  a  value  repre¬ 
sentative  of  a  He-Ne  LASER  illuminating  silicon  at  a 
temperature  of  500  K.  In  Fig.  3(a)  the  diode  is  unillumi¬ 
nated.  The  device  conductance  is  seen  to  be  negative  and 
decreasing  in  magnitude  with  increased  voltage  swing,  as 
predicted  by  analytic  theory,  equations  (16)  and  (17). 

In  Fig.  3(b)  the  diode  has  been  illuminated  from  the  n* 
region  so  that  the  majority  of  the  optically  generated  car¬ 
riers  entering  the  avalanche  zone  are  holes.  The  magni¬ 
tude  of  the  device  conductance  is  reduced,  and  this  re¬ 
duction  is  more  pronounced  at  the  higher  voltage  swings, 
in  accordance  with  (9)  and  (17).  Indeed,  for  a  peak  RF 
voltage  of  40  V,  oscillatory  solutions  are  no  longer  pos¬ 
sible:  this  is  the  basis  of  the  optical  switching  effect  [6], 
[8],  [9], 

In  Fig.  3(c)  the  direction  of  the  illumination  has  been 
reversed,  so  that  the  majority  of  the  optically  generated 
carriers  entenng  the  avalanche  zone  are  electrons.  This 
results  in  a  much  greater  admittance  change  than  that  for 
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Pig  4  W-Band  IMPATT  oscillator  equivalent  circuit  used  in  iumpuiri 
modeling  studies 

n  *  illumination  due  to  the  higher  ionization  coethcient  tor 
electrons  in  silicon  Thus  supenor  optical  control  perfor¬ 
mance  would  be  anticipated  from  this  illumination  config¬ 
uration. 

B  Optical  Tuning  Studies 

In  order  to  obtain  a  direct  comparison  with  the  optical 
tuning  theory  of  Section  Il-C  the  IMPATT  diode  model 
was  coupled  to  the  simple  oscillator  circuit  model  shown 
in  Fig  4  Lp  represents  the  diode  bond  wire  inductance. 
Rp  the  diode  parasitic  resistance,  /?,  the  load  resistance. 
C  the  dc  block  capacitance  and  Z,  a  quarter-wave  trans¬ 
former  If  Lp  is  chosen  to  resonate  with  the  capacitive  sus 
ceptance  of  the  diode  the  impedance  seen  looking  towards 
the  diode  at  plane  X  -  X’  is  real  and  negative  The  am¬ 
plitude  of  the  oscillation  is  then  determined  by  selecting 
Z,  to  match  RL  to  the  impedance  at  X  -  X '  for  the  required 
operating  point.  For  comparison  with  the  enperimental 
work  of  Section  IV  the  diode  modeled  in  this  pari  of  the 
study  was  a  W-Band  p*-n-n*  silicon  device,  with  oscil¬ 
lator  circuit  values  chosen  to  give  an  oscillation  frequency 
of  94  GHz 

To  illustrate  the  maximum  variation  in  tuning  slope  with 
optically  generated  current  composition  two  extreme  il¬ 
lumination  configurations  were  used.  In  the  first,  elec¬ 
trons  were  injected  into  the  depletion  zone  from  the  p* 
region  whereas  in  the  second  holes  were  injected  from  the 
n*  region.  Fig  5  shows  the  optical  tuning  characteristics 
obtained  from  optically  generated  current  densities  of  up 
to  2  A  cm  ,  with  a  bias  currenl  density  of  28  45  kA 
cm The  change  in  frequency  with  injected  current  is 
linear,  as  expected  from  (20).  and  the  tuning  slope  tor 
electron  injection  is  some  4  5  times  that  for  hole  injec 
tion.  Also  shown  in  Fig  5  is  the  slope  calculated  from 
(23)  This  lies  between  the  two  modeled  curves  owing  to 
the  assumption  of  equal  ionization  coefficients  in  the  an 
alytic  theory. 

C.  Transient  Study 

One  of  the  attractions  of  optical  tuning  is  that  fast  tun¬ 
ing  should  be  possible  owing  to  the  elimination  of  var¬ 
actor  charging  time  and  similar  delays  Since  (he  com¬ 
puter  model  operates  in  the  time-domain  transient  studies 
are  readily  earned  out  Fig  6  shows  the  tuning  transient 
resulting  from  a  step  change  in  optically  injected  electron 
current  from  0-2  A  cm  ' J  in  the  W-Band  oscillator  mod¬ 
eled  previously  Ninety  percent  of  the  frequency  change 
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Kig  b  Modeled  optical  tuning  transient  for  W  Band  IMPATT  oscillator 
/,.  -*  2 8  kA  cm  J#  -  0  -•  2  A  cm  ; 

was  achieved  in  55  ps.  which  compares  favorably  with 
conventional  techniques  This  may  afford  novel  systems 
possibilities,  subject  to  the  fundamental  limitation  on  tun 
ing  speed  due  to  the  oscillator  Q  factor 

IV  Experimental  Work 

Previous  work  on  the  optical  control  of  IMPATT  os 
cillators  has  mostly  been  earned  out  at  frequencies  below 
18  GHz  |8|-|  1 3)  The  object  of  this  section  is  to  describe 
recent  expenmen'al  work  on  the  optical  tuning  of  W-Band 
(  -94  GHz)  IMPATT  oscillators. 

The  diodes  used  in  this  work  were  silicon  single-dnft 
p‘-n-n*  structuns  fabricated  at  GEC  Hint  Research  Lab¬ 
oratories  These  were  mounted  in  quartz  ring  packages 
and  used  in  resonant  cap  type  oscillator  circuits,  having 
external  Q  factors  of  about  200 

The  optical  control  signal  was  generated  by  a 
GaAs/ GaAIAs  laser  emitting  at  a  wavelength  of  850  nm 
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A  50/  125-jim  graded  index  optical  tiber  was  used  to  cou¬ 
ple  the  output  to  the  edge  of  the  1MPATT  diode,  as  shown 
in  Fig  7  The  position  of  the  fiber  was  adjusted  to  give 
maximum  optically  generated  current  using  a  microma- 
mpulator  With  this  arrangement  optically  generated  re¬ 
verse  saturation  currents  of  up  to  20  5  uA  were  obtained 
for  an  optical  power  of  3  mW  The  optically  generated 
current  can  be  related  to  the  optical  power  by 


Ao 


(27) 


where  P,v  is  the  optical  power,  ij„  the  ino-mal  quantum 
efficiency,  \  the  optical  wavelength,  r  the  electronic 
charge,  r  the  velocity  of  light  in  vacuo,  and  C  a  dimen 
sionless  coupling  efficiency  Assuming  ij,  “  I  at  the  op¬ 
erating  wavelength  jives  C  -  1  pereem  This  low  value 
is  due  partly  to  refraction  in  the  quartz  ring  package  and 
partly  to  the  small  height  of  the  IMPATT  sinicture,  around 
3  urn.  compared  with  the  fiber  core  diameter  of  50  |im 
Fig  8  shows  a  typical  spectrum  analyzer  display  for  an 
oscillator  bias  current  of  ItX)  mA.  giving  a  tree-running 
frequency  of  91.83  GHz  The  center  peak  is  the  output  of 
the  umlluminated  oscillator,  while  the  right  hand  peak  is 
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the  output  ot  the  oscillator  with  an  optically  generated 
current  of  20  3  /iA.  giving  a  frequency  shift  of  9.4  MHz 
Fig  9  is  a  plot  of  oscillation  frequency  against  optically 
generated  current  The  shift  in  frequency  is  proportional 
to  injected  current  as  would  be  expected  from  (20)  with 
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small  inicclcJ  current.  The  tuning  slope  can  he  compared 
with  analytic  theory  (23k  having  estimated  the  avalanche 
gain  factor  Ftl  [  X  }  from  the  measured  output  power  ( 19}. 
For  r(J  -  0  "  ps  and  Q  -  20()  the  predicted  tuning  slope 
is  550  kHz  •  ja A ~ 1 .  From  Fig.  9  the  measured  value  is 
490  kHz  pA  1  There  is  thus  good  agreement,  particu¬ 
larly  bearing  in  mind  the  equal  ionization  coefficients  ap¬ 
proximation  inherent  in  the  theory.  These  tuning  slopes 
are  considerably  less  than  those  in  the  computer  modeling 
work  of  Section  III  B  because  of  the  low  Q  oscillator  cir¬ 
cuit  modeled.  Fig.  4. 

V.  Conclusion 

An  analytic  theory  for  the  optical  control  of  IMPATT 
diodes  has  been  described.  This  theory  has  been  supple¬ 
mented  by  results  from  large  signal  computer  modeling 
studies  Recent  experimental  measurements  on  optical 
lumng  in  W-Band  IMPATT  oscillators  have  been  pre¬ 
sented.  showing  good  agreement  with  the  theoretical 
work. 

The  computer  modeling  studies  have  shown  the  limi¬ 
tations  of  the  equal  ionization  coefficients  assumption  used 
in  the  analytic  theory.  However,  the  theory  remains  use¬ 
ful  in  identifying  the  effect  of  device  and  circuit  parame¬ 
ters  on  optical  control  performance.  The  analytic  theory 
assumes  that  the  optically  generated  circuit  is  sufficiently 
small  not  to  alter  the  diode  bias  voltage  or  oscillation  am¬ 
plitude  significantly.  Chiu  and  Freyer  [II]  have  shown 
that  where  this  condition  is  not  met  significant  increases 
m  optical  tuning  sensitivity  arise  from  thermal  effects 

The  computer  modeling  study  has  indicated  that  optical 
tuning  is  capable  of  rapid  frequency  scanning,  a  feature 
of  considerable  systems  interest  Experimental  work  is 
currently  under  way  to  test  this  prediction. 

Optical  tuning  of  a  W-Band  IMPATT  oscillator  has 
been  demonstrated  over  a  range  of  10  MHz  The  main 
cause  of  this  limited  tuning  range  is  the  poor  coupling 
efficiency  between  the  optical  source  and  the  IMPATT 
diode  An  important  area  for  further  work  is  the  fabrica¬ 
tion  of  special  millimeter-wave  IMPATT  structures,  hav¬ 
ing  an  improved  optical  response,  such  as  have  already 
been  developed  at  lower  frequencies  (21)  If  the  coupling 
efficiency  could  be  increased  to  65  percent,  a  typical  value 
for  h'gh  speed  photodiodes,  the  tuning  range  would  in¬ 
crease  to  over  600  MHz  in  the  oscillator  described. 

Optical  injection  locking,  as  described  in  Section  II -D, 
has  been  demonstrated  at  microwave  frequencies  using 
directly  modulated  semiconductor  lasers  as  optical  sources 
(12).  1 13]  However  the  maximum  modulation  frequency 
is  limited  by  earner  and  photon  lifetime  considerations 
leaser  heterodyning  (24)  offers  •  technique  for  extending 
optical  injection  locking  into  the  millimeter-wave  fre 
quency  range  if  the  formidable  technological  difficulties 
can  be  overcome 

It  is  hoped  thai  the  work  described  has  illustrated  the 
potential  of  optical  control  Although  much  further  work 


remains  to  be  done,  ihe  attractiveness  of  optical  signal 
processing  and  distribution  should  ensure  a  place  for  such 
techniques  in  future  microwave  and  millimeter-wave  s>s- 
tems. 
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